Direct measurements of the penetration depth in a superconducting film using 

magnetic force microscopy 
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We report the local measurements of the magnetic penetration depth A in a superconducting Nb 
film using magnetic force microscopy (MFM). We developed a method for quantitative extraction 
of the penetration depth from single-parameter simultaneous fits to the lateral and height profiles 
of the MFM signal, and demonstrate that the obtained value is in excellent agreement with that 
obtained from the bulk magnetization measurements. 



A fundamental property of superconductors is the abil- 
ity to expel an external magnetic field (Meissner effect), 
which is screened on the scale of the magnetic penetra- 
tion depth, A. In type II superconductors above a lower 
critical field (i?ci), however, it is energetically favorable 
to allow partial penetration of the magnetic field in quan- 
tized units called vortices. Each vortex carries one flux 
quantum, $0, and is surrounded by supercurrents, which, 
in the dilute vortex limit, decay on the same length scale, 
A. The value of A is related to the density of the su- 
perconducting electrons, and its quantitative determina- 
tion is important for understanding of superconducting 
materials, symmetry of the superconducting state and 
underlying mechanism of superconductivity. P, 0, 0, Q 
Experimental techniques that measure A in bulk samples 
include microwave measurements @], microstrip res- 
onators 01; two-coil mutual inductances Q, muon-spin- 
rotation [9| and superconductin g qu antum interference 
device (SQUID) magnetometry. [10| Recently, scanning 
probes have been app lied to in-situ measurements of A. 
Scanning SQUID [ll| and scanning Hall microscopy 12 1 
measure the vortex magnetic field above the surface of a 
sample, but the data analysis needs to take into account 
many parameters, and multistep lithography is required 
for the optimized performance and sensitivity of both 
probes, [isl. [li. [ist 



Magnetic Force Microscopy (MFM), a well established 
scanned probe technique, has been used for imaging and 
manipulation of individual vortices in thin films and sin- 
gle crystals [H, [H [13, [H, but, until now, MFM has 



not been implemented for measurements of the absolute 
values of A due to complexity of acquired images. Rose- 
man et al. [iOl estimated the penetration depth from the 
width of the lateral constant height scans of the magnetic 
force across a vortex, treating the tip as a point object, 
and found the value of A several times higher than that 
obtained by other methods. In this Letter, we report 
the measurements of the absolute value of A using the 
MFM technique. We use the knowledge of the size and 
magnetic properties of the probe tip to develop a model 
that allows us to fit the acquired MFM spectra with the 
penetration depth A as the only fitting parameter. We 



compare obtained best fitting values of A, with those mea- 
sured in a SQUID magnetometer, and demonstrate that 
our experimental and modeling approach allows for the 
extraction of the penetration depth from the MFM mea- 
surements on an individual vortex. 
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FIG. 1: (Color online) a) 1.9x1.9 fira^ scan in a Nb film at 
4.3 K. The field of view contains 53 vortices, and the vortex 
density agrees well with the expected value in 30 mT field, b) 
image of an individual vortex acquired in the external field of 1 
mT, the color scale corresponds to 0.7 Hz cantilever resonance 
frequency change which is a function of the field above the 
vortex c) cross section of the vortex and the background noise. 

All measurements described in this Letter were per- 
formed in a home-built low temperature MFM appara- 
tus at T = 4.3 K. Details of the experimental setup can 
be found in Ref. [2l|. We imaged vortices in a 300±5 
nm thick niobium film fabricated by electron beam de- 
position. With a tip-sample distance of ~ 25 /xm, the 
sample was cooled through Tc (Tc = 8.6 K measured 
with the SQUID magnetometer) to 4.3 K in the pres- 
ence of a magnetic field. The tip was then brought close 
to the sample for imaging in the dynamic MFM mode. 
We used high resolution SSS-QMFMR cantilever with a 
resonant frequency of /o ~ 70 kHz and a spring con- 
stant of c = 1.6 N/m available from Nanosensors Inc. 
(www.nanosensors.com). Fig. [T] demonstrates imaging of 
vortices in the Nb film with panel a) taken in the sample 
cooled in a field of 30 mT, and panel b) giving the pro- 
file of an individual vortex imaged in the field of 1 mT 
(yielding large intervortex spacing). It is important to 
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mention that standard MFM cantilevers have difficulties 
resolving individual vortices in fields higher than 10 mT 
17| and an ultra-sharp double-pyramid design of MFM 
tips has been used in our probes. Panel c) shows the 
cross section of the vortex, and is a typical curve used 
for fitting below. We also recorded MFM spectra as a 
function of the probe-sample separation. 
The measured frequency shift, 



A/ 



2c dz 



(1) 



depends on the force on the cantilever due to the spatially 
varying magnetic field of the vortex. At each point with a 
local magnetic moment m(r), the force F(r) — V[m(r) ■ 
B(r)], so that the net force on the tip is 



F = / V[m(r) • B(r)]dr. 

J tip 



(2) 



Previous work showed that, if the tip is treated as a point 
object, the resulting value of the penetration depth is 
nearly two times larger than expected . Consequently, 
we parameterize the tip as two cones as shown in the in- 
set of Fig. [2]d), based on scanning electron microscope 
images of a tip. The tip is covered with a magnetic Co- 
Cr film of thickness 12 nm, and the film is saturated in 
an external field prior to measurements to the satura- 
tion magnetization of Msat ~ 0.8 T [24I, which gives the 
magnetic moment per unit volume of 6.4 x 10^ A/m. 

Since the thickness of the superconducting sample is 
several times the penetration depth (the assumption con- 
firmed by the value of A wc find below), the field at dis- 
tance z above the film surface is very close to that of 
the vortex outside of a bulk superconductor, and is very 
well approximated by the field of a magnetic monopole of 
magnitude 2$o, where $0 is located at distance d below 
the surface [23|. The far, z 3> A (near, z < A) field is 
best described by = A {d ~ 1.27A) [i^, and we use the 
interpolation formula d = A(l -I- 0.27/ (1 -I- z^/A^)) in our 
fitting procedure. Importantly, with this interpolation, 
the penetration depth is the only fitting parameter. 

The fits to the MFM scans are shown in Fig. [2l 
and demonstrate that Xf ~ 109 nm gives a high qual- 
ity agreement with the experimental results. To ob- 
tain this value wc used the Lcvenbcrg-Marq uardt non- 
linear least squares fitting algorithm 2J| , where the vari- 
ance of the experimental data relative to the best fit is 
cr2 = Y,, ifi - fm{xi,Xf)f, where is the model 
prediction for a given data point (xi) and fi is the ex- 
perimental data point. Fig. [3] shows the measure of the 
fit quality, = {^(^IT'^ J2i ift ~ fm{xi, Xf)f, so that 
~ 1 means a good fit. The deeper minimum in Fig. [3] 
b) indicates that the vertical scan is more restrictive for 
our fitting procedure. Consequently, the simultaneous fit 

to both scans yields A/ = 109 ±11 nm, with the variance 

-2 

. This value exceeds the pene- 
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FIG. 2: (Color online) Comparison of the model to the ex- 
perimental data (dots) for lateral scans, a)-c), and the probe- 
sample separation (vertical) scans, d). A = 60, 109, 160 nm 
correspond to dotted, solid, and the dashed lines respectively. 
The probe-sample distance is indicated in each panel. 
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FIG. 3: Fit quality for the full width half maximum 
(FWHM) of the a) lateral scan and for the b) vertical scan 
for 26, 78, 183, and 273 nm probe-sample separations. 



tration depth for pure Nb, Aq = 39 nm, likely due to the 
presence of disorder, see below. 

We used a SQUID magnetometer to verify our result 
for A by dc magnetization. Measurements were per- 
formed on a 5x3 mm^ Nb film with the field applied nor- 
mal to the surface. A typical magnetization hysteresis 
loop M(II), obtained at T = 4.5 K after cooling the sam- 
ple from above Tc in zero field, is shown in Fig. |4l Hd 
cannot be used for a reliable determination of A since de- 
magnetization effects in this configuration are large and 
the vortex penetration field is strongly increased by pin- 
ning. In contrast, the measurement determines the upper 
critical field, Hc2 = 'I'o/27r^^(r) (marked by the vertical 
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arrow in Fig. [4]), where ^ is the superconducting coher- 
ence length. We fit i/c2(T) with a straight hne with 
= 8.6 K and the slope dH^2/dT = - 3200 Oe/K. Us- 
ing the dirty hmit results ^(T) = 0.855 v^^oV(l - T'/TJ, 
and A(r)=Ao(T) ^(^o/l-330 with the clean Nb val- 
ues ^0 = 38 nm and Aq = 39 nm, we obtain the elec- 
tronic mean free path Ik, 4.2 nm and A(0)Ri 102 nm 
(thus A(4.3K) K, 105 nm), in excellent agreement with 
the MFM fits above. 
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FIG. 4: (Color online) Measurements of the dc magnetization 
in the Nb film. Top panel: Hysteresis loop. Segment 1 is the 
Meissner response, segments 2 and 4 are in the vortex state, 
with hysteresis due to pinning. The vortex signal disappears 
along segment 3, above the upper critical field Hc2- Middle 
panel: expanded high field region used for determination of 
Hc2 at different temperatures. Bottom panel: temperature 
dependence of Hc2{T) 

To conclude, we have developed a method for reliable 
extraction of the numerical values of A from MFM images 
of vortices in superconducting films. We demonstrated 
the viability of the method by studying the Nb film. We 
expect that this approach will be used for local deter- 
mination of the penetration depth in a variety of novel 
superconductors and opens new avenues for the applica- 
tion of MFM in quantitative studies of materials. 

We acknowledge technical assistance of J. K. Baldwin 
with fabrication of Nb films. This work was supported by 
the US DOE at Los Alamos and via Grant No. DE-FG02- 
08ER46492 (I. V. and J. P. T.), by the Louisiana Board 
of Regents (J. P. T.), and was performed, in part, at the 
Center for Integrated Nanotechnologies at Los Alamos 
and Sandia National Laboratories. 



[1] R. Prozorov, R. W. Giannetta, P. Fournier and R. L. 

Greene, Phys. Rev. Lett. 85, 3700 (2000). 
[2] I. Bonalde, Brian D. Yanoff, M. B. Salamon, D. J. Van 

Harlingen, E. M. E. Chia, Z. Q. Mao and Y. Maeno, 

Phys. Rev. Lett. 85, 4775 (2000). 
[3] W. N. Hardy, D. A. Bonn, D. C. Morgan, R. Liang and 

K. Zhang, Phys. Rev. Lett. 70, 3999 (1993). 
[4] N. Schopohl and O. V. Dolgov, Rev. Lett. 80, 4761 

(1997). 

[5] W. N. Hardy, S. Kamal, and D. A. Bonn, in The Gap 
Symmetry and Fluctuations in High-Tc Superconduc- 
tors, edited by J. Bok, G. Deutscher, D. Pavuna, and 
S. A. Wolf, Plenum, New York, 1998. 

[6] C.P. Bidinosti and W.N. Hardy, Rev. Sci. lustrum. 71, 
3816 (2000). 

[7] B.W. Langley, S. M. Anlage, R. F. W. Pease, and M. R. 

Beasley, Rev. Sci. lustrum. 62, 1801 (1991). 
[8] A. T. Fiory, A. F. Hebard, P. M. Mankiewich and R. E. 

Howard, Appl. Phys. Lett. 52, 2165 (1988). 
[9] J. E. Sonier, R. F. Kiefl, J. H. Brewer, D. A. Bonn, J. F. 

Carolan, K. H. Chow, P. Dosanjh, W. N. Hardy, Ruixing 

Liang, W. A. MacFarlane, P. Mendels, G. D. Morris, T. 

M. Riseman, and J. W. Schneider, Phys. Rev. Lett. 72, 

744 (1994). 

[10] L. Civale and F. de la Cruz, Phys. Rev. B 36, 3560 
(1987), J.R. Thompson et al., Phys. Rev. B 41, 7293 

(1990) . 

[11] C. W. Hicks, T. M. Lippman, M. E. Huber, J. G. Ana- 
lytis, J-H. Chu, A. S. Erickson, I. R. Fisher and K. A. 
Moler, arXiv:0903.5260 

[12] S. J. Bending and A. Oral, J. Appl. Phys. 81, 3721 
(1997). 

[13] N. C. Koshnick, M. E. Huber, J. A. Bert, C. W. Hicks, J. 

Large, H. Edwards, and K. A. Moler, Appl. Phys. Lett. 

93, 243101 (2008). 
[14] A. M. Chang, H. D. Hallen, L. Harriott, H. F. Hess, H. 

L. Kao, J. Kwo, R. E. Miller, R. Wolfe, and J. van der 

Zielb, Appl. Phys. Lett. 61, 1974 (1992). 
[15] J. W. Guikema, H. Bluhm, D. A. Bonn, R. Liang, W. N. 

Hardy, and K. A. Moler, Phys. Rev. B 77, 104515 (2008). 
[16] A. Moser, H. J. Hug, I. Parashikov, B. Stiefel, O. Fritz, 

H. Thomas, A. Baratoff, and H.-J. Gitntherodt, Phys. 

Rev. Lett. 74, 1847 (1995). 
[17] A. Volodin, K. Temst, C. Haesendonck, Y. Bruynseraede, 

Physica B, 284 815 (2000). 
[18] M. Roseman, P. Grittter, Appl. Surf. Sci., 188 416 

(2002). 

[19] O. M. Auslaender, L. Luan, E. W. J. Straver, J. E. 
Hoffman, N. C. Koshnick, E. Zeldov, D. A. Bonn, R. 
Liang, W. N. Hardy, and K. A. Moler, Nature Phys., 5 
35 (2009). 

[20] M. Roseman, P. Grutter, New. J. Phys., 3 24 (2001). 
[21] E. Nazaretski, K. S. Graham, J. D. Thompson, J. 

A. Wright, D. V. Pelekhov, P. C. Hammel, and R. 

Movshovich, Rev. Sci. Instrum. (in press) 
[22] T. Takahashi, N. Ikeda, M. Naoe, J. Appl. Phys., 70 6056 

(1991) . 

[23] G. Carneiro and E.H. Brandt, Phys. Rev. B 61, 6370 
(2000). 

[24] W. H. Press, S. A. Teukolsky, W. T. Vetterling, and Brian 
P. Flannery, Numerical Recipes (Cambridge, 2007). 



M. Tinkham, Introduction to Superconductivity, 
Graw, Inc. (1975). 



